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We investigate the ground state persistent spin current and the pair entanglement in one- 
dimensional antiferromagnetic anisotropic Heisenberg ring with twisted boundary conditions. Solv- 
ing Bethe ansatz equations numerically, we calculate the dependence of the ground state energy 
on the total magnetic flux through the ring, and the resulting persistent current. Motivated by 
recent development of quantum entanglement theory, we study the properties of the ground state 
concurrence under the influence of the flux through the anisotropic Heisenberg ring. We also include 
an external magnetic field and discuss the properties of the persistent current and the concurrence 
in the presence of the magnetic field. 



I. INTRODUCTION 

Transport properties of strongly correlated systems 
have attracted great theoretical and experimental in- 
terest for more than two decades. In particular, low- 
dimensional systems show significant deviations in trans- 
port properties from the usual Fermi-liquid quasiparticle 
description. It has been revealed that electron correlation 
and topology play important roles in one-dimensional 
systems. Recently, the study of transport properties in 
integrable models has been an active field of research. 
The two most studied models are the hubbard model and 
the spinless fermion(or, equivalently, spin- 1/2 Heisenberg 
chain) model. Several experiments have 
observed persistent currents in mesoscopic metallic and 
semiconducting rings pierced by a magnetic flux. These 
studies have led to many theoretical investigations focus- 
ing on the interplay of electron-electron interaction and 
disorder in such systems. The persistent current in a fer- 
romagnetic Heisenberg ring has been studied in crown- 
shaped magnetic field []j , which can also be driven by an 
inhomogeneous electric fields @ due to the Aharonov- 
Casher effect |9(. Based on a spin- wave approach, the 
spin current in an antiferromagnetic Heisenberg ring with 
integer spin in an inhomogeneous magnetic field has been 
investigated very recently |l0j| . 

In the first part of this paper we study the ground state 
persistent current of a spin- 1/2 anisotropic Heisenberg 
spin ring pierced by a flux. By solving Bethe ansatz equa- 
tions numerically, we find that the increasing anisotropy 
reduces the amplitude of the persistent current. This 
phenomenon has also been noticed by G.Bouzerar [ll| 
using the Lanczos method. We also study the influence 
of an external magnetic field on the ground-state persis- 
tent current, and find that the increasing magnetic field 
reduces the amplitude of the persistent current. When 
the magnetic field is larger than a certain value (h = 3.0 
for A = 0.5, when the spins are fully polarized), the per- 
sistent current vanishes. 

Quantum entanglement, as exemplified in the singlet 
state of two spin-1/2 particles (l/v^)(| U> -| !t>)> 
is a correlation between two quantum subsystems. It 
bears some resemblance to classical correlation, but dif- 



fers in important aspects, highlighted by the observation 
of the violation of Bell's inequalities in entanglement sys- 
tems. For a special case of three binary quantum ob- 
jects(three qubits), a quantitative extension of these in- 
equalities has been proposed in terms of a measure of 
the entanglement called the "concurrence" , which takes 
values between zero and one. Several theoretical stud- 
ies 0, 0, 0] have pointed out, for example, that the 
square of the concurrence between qubits A and B, plus 
the square of the concurrence between A and C, cannot 
exceed unity [l5| . Meanwhile, it has been found exper- 
imentally that entanglement is crucial to describe mag- 
netic behavior in a quantum spin system fTfj| |. Therefore, 
we hope that the study of the entanglement in Heisenberg 
models will enhance our understanding of the quantum 
features of magnetic systems. 



We may regard a spin-1/2 chain as a collection of in- 
teracting qubits. This connection has motivated us to 
carry out the second part of the study to investigate the 
entanglement in anisotropic spin chains. There have al- 
ready been several studies on the entanglement in spin 
chains [TtI ITsl TTgt l20t l2lL I22I |23 | . In particular, recent 
efforts {24l l25l l26l [27j have been made to understand 
the quantum entanglement in the ground states of some 
many-body models, in the hope that the study of en- 
tanglement can provide a new insight into the quantum 
phase transition [2^, |2i| in these systems. Motivated by 
a recent work about the entanglement generation in per- 
sistent current qubits |30| . To this end, we study the 
properties of the ground state concurrence and its con- 
nections with persistent current. Our numerical results 
show that the concurrence is a periodic function of <j>. 
The locations of the minimum concurrence correspond 
to those of the maximum persistent current where an en- 
ergy level crossing occurs. In the presence of the external 
magnetic field, the concurrence increases with the mag- 
netic field, though the minima become less sharp as the 
magnetic field increases. 
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II. THE MODEL AND ITS SECULAR 
EQUATION 

The Hamiltonian of an anisotropic Heisenberg ring 
with twisted boundary conditions reads: 
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(1) 

where N is the number of sites. The total magnetic 
flux, or the Aharonov-Bohm flux, through the ring is $, 
and <j) — 27r$/$o where $o — hc/e is the flux quantum. 
5j , Sf , and Sf are spin- 1/2 operators at site I. By us- 
ing Wigner- Jordan transformation [3l| , the Hamiltonian 
Eq. Q can be mapped to a spinless fermion model: 

H = -i^cfCi+te-W" + C+ +1 C^ N ) 



-Vj2fa- 5)^1-5). ( 2 ) 



where J/2 can be interpreted as the hopping integral 
and V = J A as the nearest-neighbor Coulomb repul- 
sion. The spinless fermionic operators cf and Ci obey 



the anticommutation relations, and 



Ci is the lo- 



cal number operator. It is well known that the ground 
state of this model has different phases: a metallic phase 
when < A < 1 and an insulating phase when A > 1, 
where A is the anisotropy, or interaction. The former 
phase is gapless while the latter gapful. 

As in Ref. , the flux can be gauged out of the Hamil- 
tonian Eq. Q, so that solving the Schrodinger equation 
in the presence of a magnetic flux with a periodic bound- 
ary condition is equivalent to that in the absence of the 
flux but with a twisted periodic condition, namely, 

ip(xi, ■ ■ -Xi + L, ■ ■ ■ ) =exp(i4>)ip(xi,---Xi,---). (3) 

We emphasize that the characteristic length L of the cir- 
cumference is a mesoscopic scale so that inelastic scatter- 
ing does not occur. This justifies the applicability of the 
Bethe ansatz approach for the present problem, because 
the ansatz embodies nondiffractive scattering. 

From the standard quantum inverse scattering meth- 
ods (QISM) 13 1 the diagonalization of Hamiltonian 
Eq. can be transformed into solving the following 
Bethe ansatz equations: 

f sinh 2 (xj + i)\ N _^ -j-f sinh ^ (xj — x t + 2i) 



\sinh \ {xj — i) 



= ' II 



sinh Z (xj — xi — 2i) 



where 7 satisfies A = cos(7). 

Taking the logarithm of Eq. J3J, we obtain 



(4) 



M 



N9 1 ( Xj , 7) = 27r/j + 4> + J2 °2 (xj - xi, 7), (5) 



where 

The ground state is described by a symmetrical sequence 
of quantum numbers around zero. They are 

M — 1 M - 3 M - 1 
Ij - { —, ■», — 2~ }• 

The energy of the ground state can be expressed as: 



E = E - JJ2 



sin 2 (7) 



cosh(7Xj) — cos(7) ' 



(6) 



In the thermodynamics limit, it can be written in an 
integral form, 



E = E - JiVsin 2 7 



Pn(x, 4>)dx 



cosh^x) — cos 7 

00 

x exp{im[Np N (x) ~ <j)}}, (7) 

m— — 00 

where p N (x l ) = (2nli + <j>)/N and E Q = JAN/4. 

At T — 0, the equilibrium persistent current can be 
given by I{4>) — —{e/K)dE{<j))/d(j) : where E(<f>) is the 
ground-state energy of Eq. (JJ, as a function of the 
boundary condition <fi [Eq. 0]. 

III. PERSISTENT SPIN CURRENT 

To obtain the ground state of the anisotropic Heisen- 
berg model, which is a spin singlet state (N = 2M), we 
solve the Bethe ansatz equations Eq. (J5J at zero tem- 
perature and calculate the dependence of the ground- 
state energy on flux. From the Eq. (J5J, we find that 
when tf> < 7T, the quantum number sequence of the low- 
est energy state should be the same as the ground state 
when <fr — 0. But if 7r < <j) < 3n, the quantum num- 
ber sequence of the lowest energy state should be Ij = 
{- M ^,- M f 1 ,-, M f A }- Likewise, for 3vr < < 5tt, 



it should be Ij = { — 



M+3 



AI+1 
2 



A/-5 



}. In other 



words, when (j> reaches (2n + 1)tt, there is an energy 
level crossing at the lowest energy level. As the flux in- 
creased, the original ground-state energy increases while 
excited energy decreases. They are the same when <fi 
equals (2n + 1)tt. Here we solve a system with 42 spins, 
and obtain the ground state energy as in Fig. QJ. We 
can find that the ground state energy is a periodic func- 
tion with respect to the flux, with a periodicity of 27r. 
Using I((f>) — —(e/h)dE(4>)/d(f), we can easily calculate 
the persistent current as shown in Fig. (0) for different 
anisotropy. I{4>) has a sawtooth-like shape, and is also a 
periodic function of flux with a periodicity of 2tt. There 
is a discontinuity where the energy level crossing occurs. 



3 



LU 



-43.230 
-43.245 
-43.260 
-43.275 
-43.290 ■ 
-43.305 K 
-43.320 



r \ 



\ 



o.o 



1.0 



45 
40 
35 h 
30 
25 
20 
15 
10 
5 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 



FIG. 1: The ground state energy of anisotropic Heisenberg 
spin ring as a function of flux when A = 0.5403. 



FIG. 3: The magnetization M of of the lowest energy state 
as a function of external magnetic field when A = 0.5. 
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FIG. 2: The persistent current in anisotropic Heisenberg 
spin ring versus flux for different anisotropic parameters A = 
0,1/2, a/3/2,1. 
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FIG. 4: The persistent current as a function of the flux under 
the different external magnetic fields when A = 0.5. 



The existence of a discontinuity at finite A comes from 
the fact that the translational invariance is persevered 
in the presence of interaction (V = J A 7^ 0). Li and 
Ma [22] considered the current of an isotropic Heisen- 
berg ring (i.e. A = 1) using the Bethe ansatz method, 
and also found the spin current as a linear function of the 
flux. In Fig. ©, we draw the persistent current with dif- 
ferent anisotropic parameters A. The current is always 
a periodic function of the AB flux, and the increasing 
anisotropic parameter A reduces the amplitude of the 
persistent current. In the regime of < A < 1, the am- 
plitude change is weak, while is was found to decrease 
more rapidly for A > 1 [TlT |. 

In the presence of external magnetic field, the hamil- 
tonian changes into H = Hq + S z , in which Hq 
is the origin hamiltonian [Eq. without the mag- 
netic field. The magnetization S z = h(N — 2M), 
H = Hq + -|(-/V — 2M). For h = 0, the ground state 
of this system is spin singlet (N=2M), while for h ^ 0, 
the ground state is no longer spin singlet, the magnetic 
field may flip some spins. For example, when h = 0.2, the 
lowest energy state of the system with 42 spins is M = 20 



rather than M = 21. In Fig. we plot the dependence 
of the ground-state magnetization as a function of the 
external field. We can then calculate the spin current. 
In Fig. we plot the persistent current at different 
magnetic fields for anisotropy A — 0.5. We find that the 
shape of the current as a function of flux is unchanged, 
but the amplitude of the spin current decreases when the 
external magnetic field increases. And when h is small 
the amplitude decreasing velocity is slower than it when 
h is much larger. After h > 3.0, the persistent current 
decreases to zero because all the spins are polarized by 
the external magnetic field. 



IV. THE PAIR ENTANGLEMENT 

We now turn to the calculation of the ground state 
concurrence C as a function of the total magnetic flux <j>. 
Because the hamiltonian is invariant under translation, 
the entanglement between any two nearest neighboring 
sites is independent of site index. When A ^ 1, Eq. 
becomes q-deformed SU(2) algebra with A = (q+q~ 1 )/2. 
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FIG. 5: The concurrence of XXZ spin ring as a function of 
A for a system with 42 spins. 
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FIG. 6: The concurrence as a function of flux when A — > 0. 



Together with the Z 2 symmetry, we have [H, S z ] = 0, 
which simplifies the reduced density matrix pm+i) of two 
neighbor sites as 



Pi(i+i) 



/ u+ 

W! z 

z* w 2 

\ u~ 



in the standard basis | Tt>> I ti>j I lt>j an d | [[>■ 
The energy of a single pair in the system is E/N = 
Tr[pi(i + i^Hi], where Hi is the part of the Hamiltonian 
between site I and I + 1, due to the translational invari- 
ance. From the definition of entanglement, we can easily 
find that the concurrence of anisotropic Heisenberg ring 
can be calculated as l35l: 



C= ^max(0,\E gs /N ■ 



AG 



u 1(1+1) 



1), (8) 



where G^ +1 j is the two-site correlation function defined 

2 dlnZ 



by Gm+i) = Tr[eyL-p(-f3H)a lz (T[ + i z ]/Z 



where (3 = 1/T, and Z is the partition function. j^ gs 



N/3 OA 
Eg 

is the ground-state energy. In order to obtain the con- 
currence, we calculate the correlation first. From the 
definition of the two-site correlation function, we know 
that once the A-dependent eigenenergy E(A) is ob- 
tained, the correlation function is simply the first deriva- 
tive of E(A)/N with respect to A. Then we obtain the 
concurrence as a function of A as plotted in Fig. JSJ. 
This result was also noticed in the Gu et al 's work 
|3(i| . They obtained an approximative function around 
the critical point A = 1 as Cq — Ci(A — l) 2 where 
C =21n2-l~ 0.386, d = 21n2-i-^-4 ^ 0.047. 
The largest value of concurrence at A = 1 is the result of 
competition between quantum fluctuation and ordering. 
This seems to be independent of the system size: five 
qubits were considered in Ref. [33], while 1280 sites were 
considered in Ref. |36l |. 

Further more, we study the influence of the magnetic 
flux. The concurrence as a function of the flux is plotted 
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FIG. 7: The concurrence as a function of flux when A 
0.707. 



in Fig. © and in Fig. J7J for two different anisotropics 
A. We find that the ground-state concurrence is also 
a periodic function of the flux. The locations of maxi- 
mum concurrence are those of minimum energy. When 
the ground-state energy increases from the minimum to 
the maximum value, the corresponding concurrence re- 
duces from the maximum value to the minimum. The 
minima of the concurrence also correspond to the energy 
level crossings of the ground state and the first excited 
state. In the presence of the external magnetic field, the 
ground state in the same fashion as Fig.Q. We calcu- 
late the ground state concurrence as a function of flux 
at different magnetic fields, and find that when the mag- 
netic field increases, the value of the concurrence also 
increases. However, the amplitude change of the concur- 
rence decreases as the magnetic field increases as we can 
see in Fig.JSJ), after we shift these different lines along y 
axis so that the maximum coincide. When the spins are 
all fully polarized, the concurrence is independent of the 
flux. 
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FIG. 8: The concurrence as a function of flux at different 
magnetic fields when A = 0.5. The lines with !i / are 
shifted along y axis in order to make the maximal value be 
the same. 



pierced by a magnetic flux <f). By solving the BAE numer- 
ically, we obtain the ground state energy, the persistent 
current, and the concurrence, which are all periodic func- 
tions of the magnetic flux. The increasing anisotropy A 
reduces the amplitude of the persistent current. The lo- 
cations of minimum concurrence are those of maximum 
persistent current where an energy level crossing occurs. 
In the presence of an external magnetic field, we find 
that when the value of the magnetic field increased, the 
amplitude of the persistent current decreases. When the 
magnetic field is larger than a certain value (h = 3.0 for 
A = 0.5, when the spins are fully polarized), the persis- 
tent current vanishes, and the concurrence is unchanged. 
The value of concurrence also increases as the magnetic 
field increased, but the amplitude change of the concur- 
rence decreases. 



V. SUMMARY AND DISCUSSION 

In summary, we have studied the properties of the 
ground state persistent current and the entanglement de- 
fined by concurrence in the anisotropic Hciscnberg ring 
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